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Abstract 

We have presented a theoretical study of electrorotation assay based on the 
spectral representation theory. We consider unshelled and shelled spheroidal 
particles as an extension to spherical ones. From the theoretical analysis, we 
find that the coating can change the characteristic frequency at which the 
maximum rotational angular velocity occurs. The shift in the characteristic 
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frequency is attributed to a change in the dielectric properties of the bead- 
ed . coating complex with respect to those of the unshelled particles. By adjusting 
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the dielectric properties and the thickness of the coating, it is possible to 



j — , obtain good agreement between our theoretical predictions and the assay data. 
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I. INTRODUCTION 

When a suspension of colloidal particles or biological cells is exposed to an external 
electric field, the analysis of the frequency-dependent response yields valuable information 
on various processes, like the structural (Maxwell- Wagner) polarization effects [0,0. The 
polarization is characterized by a variety of the characteristic frequency-dependent changes. 
While the polarization of biological cells can be investigated by the method of dielectric spec- 
troscopy H], conventional dielectrophoresis and electrorotation (ER) analyze the frequency 
dependence of translations and rotations of single cells in an inhomogeneous and rotating 
external field, respectively fH|. With the recent advent of experimental techniques such 
as automated video analysis || as well as light scattering methods ||, the cell movements 
can be accurately monitored. In ER, the frequency dependence of the polarization leads 
to a phase shift between the induced dipole moment and the rotating field, giving rise to a 
torque acting on the particle which causes the rotation of the individual cell. 

The phenomenon of ER can be developed into a useful technique known as ER assay. The 
ER assay combines antibody technology with ER to detect analytes in aqueous solutions. 
The analyte to be detected is bound to a latex bead of known dielectric properties to form 
the analyte-bead complex, which causes a change in the dielectric properties. The change 
can be detected by ER technique, thus allowing the rapid and accurate detection of analytes 
in aqueous solutions. This method can be used to detect various analytes, the selection of 
which can be controlled by the proper choice of binding agents. 

In this work, we propose the use of the spectral representation |7j for analyzing the 
electrorotation of particles in suspensions. The spectral representation is a rigorous math- 
ematical formalism of the effective dielectric constant of a two-phase composite material 
0. It offers the advantage of the separation of material parameters (namely the dielectric 
constant and conductivity) from the cell structure information, thus simplifying the study. 
From the spectral representation, one can readily derive the dielectric dispersion spectrum, 
with the dispersion strength as well as the characteristic frequency being explicitly expressed 



in terms of the structure parameters and the materials parameters of the cell suspension (see 
section II. B below). The actual shape of the real and imaginary parts of the permittivity 
over the relaxation region can be uniquely determined by the Debye relaxation spectrum, 
parametrized by the characteristic frequencies and the dispersion strengths. So, we can study 
the impact of these parameters on the dispersion spectrum directly. The same formalism 
has been used recently to study the dielectric behavior of cell suspensions f| . 

In this connection, we mention alternative methods, namely, by solving the Laplace's 
equation directly and Gimsa's approach based on equivalent circuits |l]|§. To our knowledge, 
none of these methods separate microstructure parameters from material information. 

II. FORMALISM 

We regard a suspension as a composite system consisting of spherical or spheroidal 
particles of complex dielectric constant Z\ dispersed in a host medium of li- A uniform 
electric field Eo = Eqz is applied to the composites along the z-axis. We briefly review the 
spectral representation theory of the effective dielectric constant to establish notations. 

A. Spectral representation 

The spectral representation is a mathematical transformation of the complex effective 
dielectric constant e e . In its original form |7j], a two-phase composite material is considered, 
in which inclusions of complex dielectric constant Z\ and volume fraction p are randomly 
embedded in a host medium of l<i- The complex effective dielectric constant e e will in general 
depend on the constituent dielectric constants, the volume fraction of inclusions as well as 
the detailed microstructure of the composite materials. 

The essence of the spectral representation is to define the following transformations. If 
we denote a complex material parameter 



then the reduced effective dielectric constant 

itf(S) = l-~, (2) 

can be written as 

^VEfTT- (3) 

where n is a positive integer, i.e., n = 1, 2, ..., and _F n and s n , are the n-th microstructure 
parameters of the composite materials J7|. In Eq. ([!]), < s n < 1 is a real number, while F„ 
satisfies a sum rule J/J: 

Y,Fn=P, (4) 

n 

where p is the volume fraction of the suspended cells. 

As a result, the spectral representation is a useful theory which helps separate the ma- 
terial property from the geometric information. In what follows, we illustrate the spectral 
representation by the capacitance of simple geometry. In particular, a parallel-plate capac- 
itor is considered as an example. We will discuss two cases, namely, the series combination 
and the parallel combination. 

In the first case, if one inserts a dielectric slab of dielectric constant ei and thickness h\, 
as well as a dielectric of £2 and thickness h^ (both of the same area A), into a parallel-plate 
capacitor of total thickness h = h± + /t2, the overall capacitance C is given by 

c- 1 = C7 1 + c 2 l 

where C\ = e\A/hi, and C2 = 62^/^2, A being the area of a plate. On the other hand, we 
may define the equivalent capacitance as C = e e A/h, where e e is the effective dielectric con- 
stant. That is, we replace the composite dielectric by a homogeneous dielectric of dielectric 
constant i e . 

Let l\ = e 2 (l — 1/5), we can express C in the spectral representation, 

Ai 2 Al 2 hi/h 2 



C 



h s — \i2Jh 
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In accord with the spectral representation, one may introduce w(s) = 1 — e e /e 2 , which is in 
fact the same as w(s) = 1 — C/Cq, where Co is the capacitance when the plates are all filled 
with a dielectric material of e 2 , namely, Cq = e 2 A/h. Thus we obtain 

w s ) = TTTh- 

s — n 2 /h 

from which we find that the material parameter is separated from the geometric parameter. 
The comparison of w(s) with Eq.(|3|) yields 

F 1 = hi/h, si = h 2 /h. 

We should remark that F-± obtained herein is just equal to the volume fraction of the dielectric 
of ei, and that s\ satisfies < si < 1, as required by the spectral representation theory. 

Next we consider the parallel combination. If one inserts a material of dielectric constant 
£\ and area w\ as well as a dielectric of I2 and area u^ (both of the same thickness h) , into 
a parallel-plate capacitor of total area A = W\ + w 2 , the overall capacitance C is given by 

C = d + C 2 

where G\ = l\W\jh, and e 2 = e 2 w 2 /h. Similarly, after introducing the effective dielectric 
constant e e , we may define the overall capacitance as C = e e A/h. 
Again, in the spectral representation, let l\ = e 2 (l — 1/s), then 

C = l -2— - i2Wl 
h hs 

Writing w(s) = 1 — C '/Co, we obtain 

( ~, Wi/A 
wis) = — z — • 
s 

From this equation, the material parameter is also found to be separated from the geometric 
parameter. It is clear that F% = W\/A, i.e., the volume fraction of the dielectric of ?i, and 
si = 0. 



B. Shell-spheroidal model 

For inclusions of arbitrary shape, the spectral representation can only be solved numer- 
ically 0. However, analytic solutions can be obtained for isolated spherical and ellipsoidal 
particles. For dilute suspensions of prolate spheroidal particles, the particles can be regarded 
as noninteracting. The problem is simplified to the calculation of s n with a single particle, 
which can be solved exactly. 

In fact, the suspension of shell-spheroidal particles dispersed in a host medium is a 
three-component system. Although the spectral representation is generally valid for two- 
component composites, we have recently shown that it also applies to composites of shelled 
spheres randomly embedded in a host medium J9j . Similarly, we will show that the spectral 
representation also applies to the suspension of spheroidal particles of complex dielectric 
constant 1\ coated with a shell of l s dispersed in a host medium of e 2 , where 

I = e + a/i2nf (5) 

where / is the frequency of the applied field. In what follows, we will show that from 
the spectral representation, one can obtain the analytic expressions for the characteristic 
frequency at which the maximum electrorotation velocity occurs. The depolarization factors 
of the spheroidal particles will be described by a sum rule 

L z + 2L xy = 1 (6) 

where < L z < 1/3 and L xy are the depolarization factors along the z- and x- (or y-) axes 
of the spheroidal particle, respectively. In fact, L z = L xy = 1/3 just indicates a spherical 
particle. 

The phenomenon of electrorotation is based on the interaction between a time-varying 
electric field E and the induced dipole moment M. The dipole moment of the particle arises 
from the induced charges that accumulate at the interface of the particle. As the prolate 
spheroidal particles are easily oriented along their long axes by the rotating field, we consider 



the orientation in which the long axis lies within the field plane [fT0|j . In fact, the extension of 



our theory to deal with oblate spheroid is straightforward. In doing so, it suffices to consider 
1/3 < L z < 1 and the formlism will remain unchanged. 

The angle between M and E is denoted by 9, where 9 = u>x time and uj = 2irf is angular 
velocity of the rotating electric field. The torque acting on the particle is given by the vector 
cross product between the electric field and the dipole moment, so that only the imaginary 
part of the dipole moment contributes to the electrorotation response. In the steady state, 
the frequency-dependent rotation speed £)(/), which results from the balance between the 
torque and the viscous drag, is given by 

Q(f) = -F(e 2 , E, rj)Im{b z (cos 2 9) + 6^ (sin 2 9}), 

= -F(e 2 , E, r])Im(b z /2 + b xy /2), (7) 

where E is the strength of the applied electric field and F is a coefficient which is inversely 
proportional to the dynamic viscosity rj of the host medium. For spherical cells, F = 
—€2E 2 /2r]. In Eq.(7), Im(- ■ •) indicates the imaginary parts of (••■), and the angular brackets 
denote a time average. Since the angular velocity of the rotating field is much greater than 
the electrorotation angular velocity, i.e., uj ^> Q, the time averages are just equal to 1/2. 



For a single shelled spheroidal particle, the dipole factor b z is given by |TT 



r 1 (e s - e 2 )[e s + L z (ei - e s )} + (ei - e s )y[e s + L z (e 2 - e s )} , > 

3 (e s - ei)(e 2 - e s )yL z (l - L z ) + [e s + (ei - e s )L z ][i 2 + (e s - e 2 )L z ] 

where y is the volume ratio of the core to the whole shelled spheroid, while b xy can be 
obtained by replacing L z with L xy in Eq.(8). 

We are now in a position to represent b z and b xy in the spectral representation. Let 
l\ = £ 2 {l — 1/s), and assume x = e s /i 2 , we obtain 

b z = N.P. + J±- (9) 

s — S\ 

where N.P. denotes the nonresonant part which vanishes in the limit of unshelled 
spheroidal cells. In Eq.(9), the various quantities are given by 



Si = 


■ 
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Fl = 


—x 2 y 



«7 



(10) 
(11) 



N p _ -L z + L z y + x(-l + 2L z + y- 2L z y) + x 2 (l-L z -y + L z y) 

a 
a = 3L z - 3L 2 Z - 3L z y + 3L 2 z y + x{3 - 6L Z + 6L 2 Z + 6L z y - 6L 2 z y) 

x 2 (3L z -3L 2 z -3L z y + 3L 2 z y) 
(3 = -L z + L 2 z + L z y - L\y + x(-L 2 z - L z y + L 2 z y) 
>y = L z -L 2 z -L z y + L 2 z y + x(l - 2L Z + 2L 2 Z + 2L z y - 2L 2 z y) 

+x 2 (L z -L\- L z y + L\y) 

Note that we have assumed i to be a real number, which will be justified below. After 
substituting e = e + a/i2irf into Eq. [9] , we rewrite b z after simple manipulations 

S- Si 1 + if/fd 

with s = (1 — £1/62) _1 and t — (1 — cri/<7 2 ) _1 , where 

56 1 = Fi S_t 

/cl 



(£- si)(s- Si)' 

1 <7 2 s(t - Si) 



2-k e 2 £(s — S\ 
Similarly, we may rewrite b xy as 

F 2 , r)V 



„ = (WJ > / + _^ ) + _^_. (14) 



Therefore, we obtain 

<5e 2 = -F 2 

/c2 



(*- s 2 )(«- S 2 )' 
1 a 2 s(t - s 2 ) 



27r e 2 t(s - s 2 ) 

Note that N.P.', s 2 and F 2 are obtained by replacing L z with L xy in the expressions for 
N.P., Si and iq, respectively. 



We have thus predicted that two characteristic frequencies may appear for unshelled or 
shelled spheroidal particles. Previous theories were often limited to spherical particles, i.e., 
L z = L xy = 1/3. Therefore, only one characteristic frequency exists. We should remark that 
even though two characteristic frequencies are predicted, only one of them is dominant (see 
below) . 

III. NUMERICAL CALCULATIONS 

The model put forward in the previous section applies to various situations such as 
biological cells and polystyrene beads, etc. Here we perform numerical calculations to inves- 
tigate the characteristic frequency. Let s — 1.1, t — —0.005, and e 2 = 80eo, where eo is the 
dielectric constant of the vacuum. In Fig.l, we investigate the effect of particle shape on s\ 
and f c i (upper panels), and s 2 and / C 2 (lower panels), for z = 2 and cr 2 — 2.9 x 10~ 5 Sm~ 1 , 
where z = 1/y, and z 1 ^ 3 > 1 reflects the thickness of the shell (or coating). As is evident 
from the figure, an increase in the dielectric constant ratio x leads to a red-shift of the 
characteristic frequency. For a certain x, a small depolarization factor, i.e., the particle is 
largely deviated from spherical shape, may yield a red-shift too. 

In Fig. 2, we investigate the effect of the depolarization factor on the quantities — Im[b z ]/2 
and —Im[b z /2 + b xy /2]. It is evident that the effect of b xy on the peak is small as L z is small, 
whereas it gives a large effect as L z is large. Generally speaking, the dipole moment along 
the x-(or y-) axis strongly affects both the location and the magnitude of the peak of rotation 
speed. For spheroidal particles, only one peak is found, and the other peak predicted by the 
theory may be too small to observe. 

In order to validate our theory, here we considered the spherical particles as a limiting 
case of our model. In Fig.3, s\ and f c are plotted versus x for a spherical bead (L z = 1/3), 
for different z as cr 2 = 2.9 x 10 -5 <S'to -1 . For large x, the shell thickness has only a minor 
effect on the characteristic frequency. Moreover, a thick shell leads to a red shift (blue shift) 
of the characteristic frequency appearing when x > 1 (x < 1). Also, all the f c predicted by 
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different z are smaller (larger) than that predicted one at z = 1 (i.e., unshelled bead) for 
x > 1 (x < 1). 

At x = 1, i.e., the shell has the same dielectric constant as the host, all the characteristic 
frequencies predicted by different thickness of shell are the same. Similar conclusions can 
be obtained not only for spherical shape, but also for prolate spheroidal shape (not shown 
here) . 

We attempt to fit our theoretical predictions with experimental data, which are extracted 
from an assay |12| . In this assay, three cases were studied, all dealing with spherical particles: 
that is, unshelled beads, beads coated with an antibody with specificity for Giardia (Shell 
1), and beads coated with an antibody with specificity for Cryptosporidium (Shell 2). The 
bead diameter is 6/im according to Burt et al [fL^] , while Giardia and Cryptosporidium 
are both 0.8 fim in diameter. Let <j 2 = 2.18 x 10~ 5 Sm -1 , z = 6, x = 7.06 (Shell 1) and 
4.63 (Shell 2), and F = 0.353 (Without shell), 1.629 (Shell 1) and 2.278(Shell 2). Good 
agreement between our theoretical predictions and the assay data is shown in Fig. 4. From 
our theory, it is easy to find the corresponding characteristic frequencies at which maximum 
rotational angular velocity occurs, f c = 4.755 x 10 5 Hz (Without shell), 10 5 Hz (Shell 1) and 
1.415 x 10 5 Hz (Shell 2). From this, we find that the coating leads to a red shift of the 
characteristic frequency. It is because the dielectric properties of the bead-coating complex 
have been changed. 

We used x > 1 in our fitting, i.e., x = 7.06 and 4.63 for bead coated with an antibody with 
specificity for Giardia and Cryptosporidium, respectively. It is known that, for biological 
cells, such as Giardia and Cryptosporidium, they have structures such as the cell wall, 
plasma membrane and the cytoplasm, among which the cell wall has a larger conductivity 
than the suspending medium. Hence, we may safely take x > 1. 
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IV. DISCUSSION AND CONCLUSION 

Here we would like to add some comments. We would like to clarify the assumptions 
with our model in more detail. In fact, there is only one (rather than two) peak for each 
polarization in our theory. It is because we assumed the ratio of the shell to host dielectric 
constant a; to be a real and positive number. If we had retained the (indeed small) imaginary 
part of x in our calculation, then we would have two peaks for each polarization. The 
conductivity-dominated peak would have occured at substantially lower frequency. Thus, the 
neglect of the imaginary part of x is to drop the lower frequency peak. Moreover, according 
to our calculations, there is one dominant pole associated with b z and two (degenerate) sub- 
dominant poles associated with b xy in the spectral representation. Thus, for shelled particles 
in the present work, only one peak has been shown. 

We developed simple equations to describe the electrorotation of particles in a suspen- 
sion from the spectral representation. These equations serve as a basis which describe the 
parameter dependence of the polarization and thereby enhances the applicability of vari- 
ous cell models for the analysis of the polarization mechanisms. In this connection, the 
shell-spheroidal cell model may readily be extended to multi-shell cell model. However, we 
believe that the multi-shell nature of the cell may have a minor effect on the electrorotation 
spectrum. 

We have considered the isolated cell case, which is a valid assumption for low concen- 
tration of cells. However, for a higher concentration of cells, we should consider the mutual 
interaction between cells. For a randomly dispersed cells suspension, we may replace the 
dielectric constant of the host medium by the effective dielectric constant of the whole sus- 
pension. 

When a strong rotating electric field is applied to a suspension, the induced dipole 
moment will induce an overall attractive force between the polarized cells, leading to rapid 
formation of sheet-like structures in the plane of the rotating field. In reality there is a phase 
shift between the induced dipole moment of the structure and the applied field, and this can 
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lead to electrorotation. However, the situation will be much more difficult than the single- 
cell case that we have studied because the many-body as well as multipolar interactions 
between the particles will produce a complicated electrorotation spectrum. Fortunately, 
our recently developed integral equation formalism jHJ can definitely help to solve for the 
Maxwell- Wagner relaxation spectrum. 

In summary, we have presented a theoretical study of electrorotation assay based on the 
spectral representation theory. We consider unshelled and shelled spheroidal particles as 
an extension to spherical ones. From the theoretical analysis, we find that the coating can 
change the characteristic frequency at which maximum rotational angular velocity occurs. 
By adjusting the dielectric properties and the thickness of the coating, it is possible to obtain 
good agreement between our theoretical predictions and the assay data. 
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FIGURES 
FIG. 1. Upper panels: si and f c \ plotted against x for different L z at z = 2. Lower panels: s2 

and f C 2 plotted against x for different L z at z = 2. 

FIG. 2. — Im(6^/2 + b xy /2) and — Imb z /2 are plotted against frequency uj for different L z at 
z = 6, x = 2 and (j 2 = 2.9 x 10 _5 S , m~ 1 . 

FIG. 3. si and f c are plotted versus x for different z as L z = 1/3 (i.e., spherical shape). 

FIG. 4. Curve fitting for L z = 1/3 (i.e., spherical particles). 



14 



<f> 









I 

L z =0.05 


n 
\\ 


L z =0.2 


\\ 

\\ 
\ \ 


L z =l/3 


; \ \ 
\ \ \ 

': \ \. 

': \ ^^ 

\ 
\ 


- 




^~~---____ 


i 







1 

x 




<r> 







\ 


L z =0.05 


\\\ 


- L z =0.2 


\\\ 

\\\ 


L z =l/3 


\\ \ 




\ \ ■■-. 




\ x 




\ X 


"•--..._ 




i 







1 

x 




Fig.l 



15 



0.5 



0.4 



0.3 



0.2 



0.1 







-lm[b/2] as L =0.05 
-lm[b/2+b Iy /2] as L I= 0.05 
a -lm[b/2] as L 2 =0.2 
-Imfb/ifb^/J] as L 2 =0.2 
-Imib^] as L =1/3 
-lm[b/2+b Iy /2] as L =1/3 




7 



Log f (Hz) 



Fig.2 



16 



c/> 



"A 



•\ w 



z=1 

z=1.5 

z=2 

z=6 

z=27 



X 



^ 



\ 



\ w 



— z=1 
z=1.5 
z=2 
z=6 
z=27 



1 
x 



1 

X 



Fig.3 



17 



<f) 





I 
Without shell (assay) 


1 ' 








D Shell 1 (assay) 










AShell2(assay) 










Without shell (theory) 










Shell 1 (theory) , ' 

Shell 2 (theory) A 


'\/ \ 












0.2 


/ 

/ 


1 \ \ 




- 




A 


i \ 


V 






I 


/ \ 








/ 


/ \ 








/ G~ / 






0.1 


/ /D / 

a// / 


\ 

\ 


°\ 


- 




// 


\ 


v A 






_ 'u' 




\ 
\ 




r\ r\ 


1 


1 


\ A 


o \ 

^ A 



4.0 



4.5 



5.0 



5.5 
Log f (Hz) 



6.0 



6.5 



7.0 



Fig.4 



